Takahashi Integral Equation 
and High- Temperature Expansion of the Heisenberg Chain 



Masahiro Shiroishi and Minoru Takahashi 
Institute for Solid State Physics, University of Tokyo, 
Kashiwanoha 5-1-5, Kashiwa, Chiba, 277-8581 Japan 
(Dated:) 

Recently a new integral equation describing the thermodynamics of the ID Heisenberg model 
was discovered by Takahashi. Using the integral equation we have succeeded in obtaining the high 
temperature expansion of the specific heat and the magnetic susceptibility up to 0(( J/T)^*^"). This 
is much higher than those obtained so far by the standard methods such as the linked-cluster 
algorithm. Our results will be useful to examine various approximation methods to extrapolate the 
high temperature expansion to the low temperature region. 
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The spin- 1/2 Heisenberg chain is the one of the fun- 
damental models which has been continuously investi- 
gated in the realm of the magnetism in low dimen- 
sions. In 1931, Bethe constructed the eigenstates for 
the isotropic case which was later generalized to 
anisotropic cases such as the XX Z and the XYZ chain 
(see 1^). The method, which is nowadays called Bethe 
ansatz method, has provided us vast information on the 
model. Especially the thermodynamics has been studied 
through the Thermodynamic Bethe Ansatz (TBA) equa- 
tions g, ||, ^ ||, 0, ItII and the Quantum Transfer Matrix 
method (QTM) 



0, |T|,[ii,[illiE10- 



Quite recently, one of the authors (MT) has found yet 
another integral equation JTzt , which determines the free- 
energy of the XX Z chain. The integral equation consists 
of one unknown function and is very simple. We have 
found that it actually gives the same numerical values 
for physical quantities as those calculated by the TBA 
equations and the QTM method. In fact the equation 
was originally discovered in an attempt to simplify the 
TBA equations ||l^. Soon after that it was derived also 
from the fusion relation of the QTM p8|. 



In this letter, as one of the applications of this inte- 
gral equation, we study the high temperature expansion 
(HTE). As is well known, there is a systematic way to 
calculate the HTE for any lattice models based on the 
linked-cluster algorithm |l^, ^ . Particularly, in the case 
of the XXX chain, the HTE has been achieved up to 
0{{J/Tf^) (20). We have succeeded in obtaining the 
HTE to much higher order, namely, up to 0((J/r)i''°) 
by use of the new integral equation. Below we shall re- 
port the result for the isotropic XXX model. The gen- 
eralization to the anisotropic models will be published in 
a separate paper. 



The Hamiltonian of the spin-1/2 Heisenberg XXX 



chain is defined by 

N r 



h=-jy: 

N 



QX QX I QV QV I QZ QZ 
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where Sj''^'^ are the local spin-1/2 operators acting on 
the site j.We assume the periodic boundary conditions 
Sn-\-i = Si . Note that in our definition, the coupling con- 
stant J is positive for the ferromagnetic case and negative 
for the antiferromagnetic case. 

Takahashi's integral equation for the isotropic XXX 
case is given by 

u{x) = 2 cosh /i/T 




X — y + 2i 
f = -Tln«(0), 



where the contour C is a loop surrounding the origin in 
a counterclockwise manner. The equation can be solved 
numerically. 

In the following we derive the HTE of u{x). Actually 
what we have to do is only to assume u{x) in the form 



u{x) — exp 



Y,a„ix) {J/TY 



= e'^«(-)|l + ai(x)J/T+ (^a2{x) + \i{xf'^ [J/Tf 
+ (az{x) + a2[x)ai{x) + ^ai{xf \ {J/Tf + ■■■ I. 



(3) 
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and substitute it into the equation . Then by compar- 
ing the same order of J /T in the LHS and the RHS, we 
can get the equations, which, to our surprise, determine 
the functions a„(x) recursively. For example, from the 
0-th order, we get immediately 



ao(x) = In (2 cosh ft. /T) . 



(4) 



Similarly by comparing the 1-st order, we have an equa- 
tion. 



ai(x) 



4cosh^ h/T [Jc \ y{y + 2i){x-y- 2i) 

-2 1 dy 

j/(y-2i)(2;-y + 2i) j 2m 



(5) 



Noting that the second term in the RHS vanishes because 
the integrand is regular at y = 0, we can calculate ai{x) 
explicitly as 



cosh^ h/Tx^ +A 



(6) 



Repeating the similar procedures we can derive each 
an{x) successively. For example, we have found 



02 (x) 



az(x) 



12 



4cosh2 h/T (a;2+4)2 



1 



a;2 -h 6 



Acosh^ h/T (x2 -h4)2' 

1 ix^ + 36a;2 + 160 
24 cosh^ h/T (a;2 + 4)3 



(7) 



1 



lla;-^ + 36 



4cosh''/i/r (a;2+4)3 
1 3a;4 + 30a;2 



80 



24cosh^ft/T (x2-H4)3 



(8) 



In fact, with the help of Mathematical we have calculated 
a„(x) up to n = 100, which provides the HTE for the 
free energy f/T = - j:Zo dnmJ/TT up to (J/r)i"o. 



Some of the lower terms are given as 
f/T - - ln(2cosh(Vr)) + ^(1 - B^) 
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3072r4 
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5120r' 



(1 - B^) 

(1 -S2)(3 + 7S2 + 10B^) 
(1 - B^){3 -B^ - 9B^ - 21B^) 
(1-B2)(1 + 2B2) 



7J6 



X (15-5^ + 215^-635^) 
.(1-5^) 



122880T6' 

X (3 - 355^ - SbB^ - 955^ - 305^ + 3305^°) 

+ ..., (9) 

where we have set B — tanh(ft,/T). From the expansion 
(^, one can get the HTE for other physical quantities. 
For example, we list coefficients of the HTE for the 

specific heat C = —T-q?^ and the magnetic susceptibility 

X — ~ fi^ at zero magnetic field in Table | and Table ||. 
Unfortunately due to the lack of space, we can present 
the coefficients only up to 0((J/T)^^). The coefficients 
of higer order will be sent on demand to any interested 
reader. We remark that the coefficients up to 0(( J/T)^"') 
completely coincide with those given in [^.( Note the 
differences of the conventions, J ^ — J and h 2h. ) 
The terms with the order larger than 0{{J/T)'^^) are our 
new results. It will probably be impossible to get the 
HTE to such a high order using a conventional method. 




FIG. 1: Specific heat for the antiferromagnetic XXX chain 
at /i = 0. 

Usually in order to extrapolate the high temperature 
expansion series to the low temperature region, some fur- 
ther approximation methods are used. Actually the orig- 
inal series will not converge in the low temperature region 
(T about less than 0.55 ), because of the existence of the 
singularities on the complex plane with respect to the 
inverse temperature. 

Here we have applied the standard Fade approxima- 
tion to our HTE. The results are shown in Fig. l,...,Fig. 4. 
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FIG. 2: Magnetic susceptibility for the antiferromagnetic 
XXX chain at /i = 0. x(0) is taken from 
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FIG. 3: Specific heat for the ferromagnetic XXX chain at 
ft = 0. 



Since we have obtained the HTE up to 0({J/Tf^^), the 
expressions up to [50, 50] Padc approximant are available. 
For comparison, we have also plotted the numerical data 
calculated by the QTM method |jll|, (Note that the 
physical quantities in Fig. l,...,Fig. 4 were first calculated 
by the TBA equations. See particularly ^ for the fer- 
romagnetic case.) From Fig. l,...,Fig. 4, we find that the 
[50, 50] Fade approximant show good coincidence with 
data by the QTM to very low temperature region some- 
what like T ^ 0.05 except for the the magnetic suscep- 




FIG. 4; Magnetic susceptibility for the ferromagnetic XXX 
chain at /i = 0. Note that the plain series of HTE deviates 
from the QTM data around T = 0.52. 



tibility for the antiferromagnetic case (Fig. 2.). In that 
case the logarithmic anomaly around T — > is so strong 
[ p^ , [Tsl , and it probably prevents the good conver- 
gence of the Fade approximation [|o|. 

Apart from the very low temperature region, we have 
found our higher order Fade approximation gives the 
physical quantities with extremely high precision. For 
example, we have estimated the peak position of the spe- 
cific heat and the magnetic susceptibility for the antifer- 
romagnetic case by use of our [50, 50] Fade approxima- 
tion. The result for the specific heat is 



^max ^ 0.3497121234553176, 
jnmax /, ji ^ 0.4802848685890477, 



and that for magnetic susceptibility is 



(10) 



^11 



jJl = 0.6408510308513831. 



(11) 



These values are perfectly identical to the ones in [ p^ , 
where the non-linear integral equations for the QTM were 
solved numerically very carefully. (Note that our 
is four times larger than that in [Q , because of different 
normalization factors.) Note also that the peak position 
of the magnetic susceptibility (^l|) was first determined 
in pi. 



For comparison we have also estimated the peak posi- 
tion of the specific heat for the ferromagnetic case as 



(^max ^ 0.1342441913136996, 
j.max^ J ^ 0.3326119630964252. 



(12) 



In conclusion, we have shown that the the new integral 
equation by Takahashi is very useful to calculate analyti- 
cally the HTE for the ID Heisenberg model. As far as we 
know, the HTE of such a high order as 0((r/ J)i°°) has 
not been achieved for any other models except for free 
ones. We have seen that our HTE data together with 
Fade approximation, provide the very accurate numeri- 
cal data for physical quantities to sufficiently low tem- 
peratures. There are several further methods to improve 
the Fade approximation, for example, to allow for the 
asyrnptotic expression as T — > 0, etc. (see, for example, 
pO| , p2[). Based on the present results we shall investi- 
gate such possibilities in the forthcoming publication. 

The authors are thankful to A. Kliimper, Y. Nishiyama 
and K. Sakai for valuable discussions. This work is sup- 
ported by Grants-in-Aid for the Scientific Research (B) 
No. 11440103 from the Ministry of Education, Culture, 
Sports, Science and Technology, Japan. 
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TABLE I: Series coefficients q;„ for the high temperature expansion of the specific heat C = ^(^)"at/i. = 0. 



n 






n 
26 


484455914465376683487755420408217600 


1 





27 


-1592964364128699671723658807556964352 


2 


6 


28 


-1659222341377723674454893065936371187712 


3 


-36 


29 


53827694891210973745020673240061454581760 


4 


-360 


30 


5090517962961447184851808942927438864711680 


5 


7200 


31 


-388446833192725659973817494776649147157053440 


6 


15120 


32 


-11028658525378274359384407389662010654796546048 


7 


-1848672 


33 


2255854109806569120380670028755308714386167693312 


8 


11426688 


34 


-18878702580622989070078793482363993425701807063040 


9 


594846720 


35 


-11721570087037734701356860480896473609272542720163840 


10 


-11558004480 


36 


527183038642469386328859769728390518893185382125404160 


11 


-199812856320 


37 


52548749252010967993948480669499712309299856992951599104 


12 


10106191180800 


38 


-5382365237582398925074954773487741035075672601159589167104 


13 


19376365252608 


39 


-150281021589219619860159284209265140804955107276364175114240 


14 


-9289795522775040 


40 


44482678475307391762958213932359681737961800852665438128046080 


15 


121944211136778240 


41 


-670778300712303276022754187872671936481343744506621812675706880 


16 


8791781390116945920 


42 


-323311185126253530334911142992092649497388429937499549362387156992 


17 


-310402124957945954304 


43 


19271391500067613736198673193545354611765664770995927250862568636416 


18 


-7225535925744106143744 


44 


1963797073102024140530884388201619017857423297479642613447074848440320 


19 


643407197363813620776960 


45 


-261757449501391383349154989821467694962901907072496780634929173522022400 


20 


96147483542540314214400 


46 


-6715036186134671522475926929150627328836680429076585020863949295740518400 


21 


-1279121513829538179364945920 


47 


2897640509780835688484069216581936412870887902144153768250804439297575354368 


22 


27962069861743501862336200704 


48 


-67884583842448252705729493380589916284543590592089243545625816663055684075520 


23 


2398518627113966015427501883392 


49 


-278396673545453495106462073222678393434494076172191005087521563009127066632192 00 


24 


-129834725539335848980192847462400 


50 


2130333568970965233678580974509426707048585535358286474483865352948915543579033600 


25 


-3493877000064415911285457158144000 


51 


216827879657653769500650387534438914339017251205042192428944 424756474567924803174400 



TABLE IL Series coefficients jSn for the high temperature expansion of the magnetic susceptibiUty x = y X^n (n+iy. 
h = 0. 



n 




n 


/3„ 





1 


26 


108250895627317866042581831619969024 


1 


4 


27 


-17922781862082598948422103131404369920 


2 





28 


120323704775766862241382488276073447424 


3 


-64 


29 


63726499196979511634418116940394187980800 


4 


400 


30 


-2379574083879736929901931437077028466065408 


5 


4032 


31 


-199921736549129208105099470038165104315334656 


6 


-89376 


32 


17205669691857844357030111005301702510893858816 


7 


-163840 


33 


408607241970908762765181373729511748361273737216 


8 


26313984 


34 


-1024()5082506734431652137262962322034251260126822400 


9 


-191334400 


35 


1315775281576476974395821022113890916515475892469760 


10 


-9565698560 


36 


545540112377157932817560176143490374271324160383778816 


11 


210597986304 


37 


-27991080178083294022121824180738137056773157996564840448 


12 


3486950684672 


38 


-2463499970265937469283180992880401729855043417164218368000 


13 


-203634731188224 


39 


284868319106208295442521316317947214008165828662443533926400 


14 


-127324657152000 


40 


6335450833027827824904755798342667046011674072210769010753536 


15 


205019990184689664 


41 


-2393256391170887534773811023482335548575993905771247989063417856 


16 


-3169755454477500416 


42 


49322825714830688230041284007303368927055837926052566121716908032 


17 


-208763541 109969256448 


43 


17655809163462320179643661041099923582156280899288659229470449729536 


18 


8342101010835559022592 


44 


-1192343745774700569487973174622411288186826884326738448519272595456000 


19 


175912858271144581529600 


45 


-106665300086792728016396630723577790891277245485607114802374328230871040 


20 


-18366266410738921187573760 


46 


16097399250631032968718675605919227751757635356066434920608362971107688448 


21 


40780317289246872850923520 


47 


306091410138596910653186173722480725145635803390769178021236250893761904640 


22 


38668138493195891009425244160 


48 


-180022014647714649285609157196673479185707562483390012375002362479536 725032960 


23 


-983734184997038611238624428032 


49 


5408449632892687473365326165267541130523650154708192975684794827854640578560000 


24 


-75650797544886562610211717120000 


50 


17426550003663298365992833617318506452913828517681 69407986516632569765546101309440 


25 


4622511990582180728868776781545472 


51 


1503021976366022386779402749953214978007681637342006785713651741321010 78803395117056 
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